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Abstract — A modular method was suggested before to recover 
a band limited signal from the sample and hold and linearly 
interpolated (or, in general, an nth-order-hold) version of the 
regular samples. In this paper a novel approach for compensating 
the distortion of any interpolation based on modular method has 
been proposed. In this method the performance of the modular 
method is optimized by adding only some simply calculated 
coefficients. This approach causes drastic improvement in terms 
of signal-to-noise ratios with fewer modules compared to the 
classical modular method. Simulation results clearly confirm 
the improvement of the proposed method and also its superior 
robustness against additive noise. 

Index Terms — Compensating distortion, interpolation, modular 
method, optimum coefficients. 



I. Introduction 

DIGITAL to analog converters are common in digital 
signal processing and communication systems to recon- 
struct an analog signal from its discrete time samples. Several 
methods with different names were introduced in the literature 
in 1970's and 1980's |1|. S&H and LI were the dominant 
methods before that time; today, Polynomial interpolation and 
B-Spline are the usual interpolation functions Q-Q. 

These interpolators create some distortion at the Nyquist 
rate after low pass filtering, especially when S&H or LI are 
utilized. The advantage of these types of interpolators is their 
simplicity which makes them proper for practical use. To 
alleviate this problem, several methods such as inverse Sine 
filtering, over-sampling, nonlinear and adaptive algorithms 
|5|-[7), a modular method of the recovery of a signal from 
its sampled-and-held version are described in |8) for the 
uniform samples, |9j for the nonuniform samples, and succes- 
sive approximation using an iterative method [:T0)-fT2) were 
introduced. The modular method is compared to the inverse 
Sine filtering in [ 8 1 which shows that by using a few numbers 
of modules, the performance of the modular method excels the 
inverse filtering as far as noise is concerned. Over-sampling 
is not a practical solution due to its bandwidth requirements. 
The iterative method [11] outperforms the modular method at 
the cost of more computation. 

We propose an optimized modular method which enhances 
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the performance of the classical modular method fl3) . p] Our 
method is based on some optimum coefficients which are 
computed very simply by solving a least square problem. 
Indeed, these coefficients are calculated just one time for a 
specific interpolation system and are independent of signal 
to which the modular method is going to be applied. The 
coefficients themselves do not increase the complexity of the 
modular method and is very simple for practical usages. The 
simulation results show that the coefficients are well optimized 
and perform better then classical method. 

The rest of this paper is organized as follows: Section [II] 
describes our general framework and introduces the terms and 



concepts used throughout the paper. Section III introduces our 
proposed method and is a straight forward manner to find 
the optimum coefficients for the modular method. Simulation 
results and comparison with the classical modular method for 
various interpolation systems will be presented in section |IV] 
and finally, section [V] will conclude this paper. 

II. Preliminaries 

In this section we give a brief overview of the modular 
method [8| that compensates the distortion of any interpolator 
such as Sample and Hold (S&H) and linear order hold by 
mixing the sum of cosine waves and then passing them through 
a lowpass filter. Suppose x(t) is sampled at the Nyquist rate 
(£) and assume s(t) is any interpolating function that fits the 
samples of x(t). According to these assumptions it can be 
formulated as follows: 



s(f) = h(f) * V x(nT)5(t - nT) 



(1) 



where h(t) is the impulse response of the interpolation func- 
tion. The above equation can be written as shown below in 
frequency domain: 



S(f)=H(f)xY j X(f-i/T) 



(2) 



According to the modular method, an improved reconstruc- 
tion of x(t) can be derived from s{t) by following process: 



.•;-(/)= Mn\l+2cos( 2 ^) + ---+2cos( 2 ^)\*n(i) (3) 



'Part of this paper is presented in 18th International Conference on 
Telecommunication (ICT) 2011. 
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Fig. 1. The reconstruction block diagram using Modular Method 



where 11(f) is a lowpass filter with a bandwidth equal to 
the bandwidth of x(t) (W). Eq. Q can be rewritten in the 
frequency domain as follows: 



Mf) = Tl(f) J] S(f-j/T) 

j=-N 

+N +oo 

= n(f) ^ H(f - j/T) 2 *(/ - «yr - y/r) (4) 



Since is band-limited Q can be simplified as follows: 



+.<v 



X(f) = x 11(f) J] //(/ - j/D 



(5) 



/=-JV 



Therefore, it is obvious that as N increases x will per- 
fectly converge to x if Il(/) Yj]=-n H(f ~ j/T) becomes unity 
which always occurs for all interpolation functions since 
r'lZjr-ooW - j/T)) = <W, for instance this summation 
becomes summation of sine functions which is unity for 
Sample and Hold (S&H) interpolation function. However, 
practically it is not possible to apply infinite numbers of 
modules and only limited numbers of them are implemented 
by means of oscillators shown in Fig. [T] So the method will 
have distortion and to measure the distortion, let us define the 



mean-square error as: 



X+W ( +N 
V H(f-j/T)-l 



df 



(6) 



and it is obvious that lim e = 0. 

N-too 

Another issue with this approach is that the signals are 
considered to be analog. In practice, most of the signals that 
we are dealing with are discrete especially in a computer for 
interpolation of images and audio files. 

In the next section, we will drive all these relations again 
in the discrete domain and also minimize the mean square 
error using the optimization coefficients. The second part is 
the main part of this paper. 



III. Proposed Optimized Modular Method 

A. Modular method in Diecrete Domain 

Consider x[n] is band-limited discrete signal which is down 
sampled and interpolated at the Nyquist rate (1/T) by h[n] 
and the result is s[n]. Therefore: 



s[n] = h[n] * x[kT]6[n - kT] 
and equivalently: 



k=— oo 



S(k) = H(k)* |] X((k-iN/T)h 



(7) 



(8) 



where S (k), H(k) and X(k) respectively are N-point DFTs of 
s[ri], h[n] and x[n]. Supposing that is divisible by T, the 
modular method can be formlated as follows: 



x[n] = LPF{s[«] (1 + 2cos{2jm/T) + ■■■ + 2cos{2Mnn/T)) } 

(9) 

where LPF is a FFT lowpass filter. Also it can be shown that 
applying more than [T/2] modules not only does not enhance 
the performance but also may distort it. Consider T is even 
integer number and we have applied k modules more than T/2, 
then we have: 



l+S£*2o«(^) 



= 1 + 5£j 2cos(^) + X ; ' ;., 2cos( 2 -p) 
= 1 + £Jj 2cos(^f) + Z k j=l 2cos{ 2 ^p^; 



1 + 2cos(^) + Eti (-l)"2««(^f ) 



linn \ 



(10) 



It is obvious that the third term in the above final result will 
distort previous modules; this effect can be shown for odd Ts 
in the same way too. Hence, the maximum number of modules 
that is able to be applied is [T/2]. 

Furthermore, in the case that the maximum number of 
modules are applied, if we put the multiplicand of last cosine 
1 instead of 2, we will reach the impulse train. By means of 
Fourier series it would be proved as follows: 



JlnknlT 



S[n] = X ^ ~ ' T) = Z 7 

x T 
772-1 

1 +(-1)" + Yj 2cos(2nin/T) 



772-1 



1+ ^ 2cos(2mn/T) + cos(2n—n/T) 



i=l 



(11) 



So, in this situation it will gather the original samples and if 
the filter is ideal, the output would be perfectly interpolated]^] 
The equation (|9]l will be rewritten in the frequency domain 
like below: 

2 Also true for analog D/A where 2, 2, . . . , 2, 1 is not optimum but better 
than before. 
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+M 

*(*) = X(k) x U(k) J] H((k - jN/T)) N (12) 

j=-M 

Simultaneously, the error of the interpolation process would 
be formulated in the same manner performed in the pervious 
section: 

+N/T I +M \ 2 

e=Y, Y, mk - jNIT))N ~ l (13) 

-N/T \j=-M 

Again minimizing this error is our main goal and is related 
directly to the performance of our system. In next subsection, 
our proposed Technique will be introduced to achieve this 
purpose. Again minimizing this error is our main goal and 
is related directly to the performance of our system. In next 
subsection, our proposed Technique will be introduced to 
achieve this purpose. 

B. Optimization Coefficients 



Our goal is to minimize (13 1 and as a result reduce 
distortion in the modular method in order to reach more 
precise interpolation for x. Our idea is that modules could be 
applied with some coefficients shown in Fig. [2] By choosing 
these coefficients appropriately, the performance of the method 
increases efficiently and we can achieve the same result with 
fewer modules. 

Now consider modules are multiplied by cf. 



2nn IMnn 
cos{—) + ■■■ + 2c M cos( - ) 



x[n] = LPf|s[«] ^1 + 2ci 
and therefore, ( fTJfr becomes: 

+N/T ( +M 



(14) 



(15) 



-NIT \j=-M 

For convenience the following parameters are defined: 



6 = 



ffi(0) H 2 (0) 
ffi(l) H 2 (l) 

Hi(N/T) H 2 (N/T) 



H M (0) ^ 

n M (i) 
H M (N/n 



(17) 







T 




( H (0) ) 


C2 




1 


1 


H (l) 




,93 = 












~ 2 




\ C M, 




X 







(18) 



Now from these new definitions the error function can be 
simply rewritten as follows: 

e = (norm(§(£ - 23)) 2 (19) 
Hence, to minimize this error, we should solve: 

£G = 93 (20) 

Since the number of equations (N/T) are much more than 
the unknowns (M), we could find the optimum answer by 
considering the mean square error method and the pseudo- 
inverse: 



(S r S) _1 S r 8 



(21) 



Therefore, the problem can be solved and we have found 
some coefficients which minimize the error for the finite 
number of modules without any exception on interpolation 
function. Furthermore, we will show in next section that 
these coefficients cause dramatic results in comparison to 
the classical modular method which is a special case of our 
method by assigning the coefficients one. The main key in our 
method is that the coefficients are calculated very easy and 
fast and by only accessing impulse response of interpolation 
function (}i[n\). Moreover the coefficients are calculated just 
one time for an interpolation function and are stored in a 
lookup table and does not need to find them again every time 
we need them. 



Hj(k) = H((k-iN/T)) N + H((k + iN/T)) N 




Fig. 2. The reconstruction block diagram using Discrete Modular Method 
and Optimazation coefficients 



C. Optimized Modular Method for 2-D Signals 

Due to increasing need to analyze the two-dimensional sig- 
nals like image and video, we want to generalize our proposed 
method to 2-D space. For this reason, consider x{t\,t 2 ) is 
band-limited discrete signal which is sampled at the Nyquist 
rate (^-, jr) Now consider that this signal is reconstructed 
by interpolating function h(t \ , t 2 ) and the resulted signal is 
s(t\,t 2 ). With no doubt the amount of signal s{t\,t 2 ) is equal 
to x(fi,?2) in sampled points. By knowing these assumptions 
we can formulate them as follows: 

s(ti,t 2 ) = h(t u t 2 ) * V x(n\Ti,n 2 T 2 )5{t\ - n\T\,t 2 - n 2 T 2 ) 

(22) 

Easily the above equation, by a simple Fourier transform 
and by means of its properties in frequency domain, can be 
written like below: 
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X(fi-h/Tuf2-ii/T 2 ) (23) e= f f f H^-hlTuh-hlTi) 



By considering the modular method for one-dimensional 
signals we can understand that the main concept in modular 
method is rooted from averaging in frequency domain. Using 
this idea in one-dimensional domain we can generalize the 
modular method to two-dimensional domain. By knowing this, 
2 - D modular method should include the shifted version of 
that in different directions and here we use Lattice pattern. 
With these information, the reconstructed compensated signal 
x(?i , t 2 ) by means of 2 - D modular method can be written 
from signal s{f\,t%) is: 



2nt\ 2nt 2 
x(ti,t 2 ) = s{x\,x 2 ){\ + 2cos{^—) + 2cos{—-) 

+4cos(^P-) + Acosi 7 ^) + ...)* 11(f) (24) 
Similar to one-dimensional signal we can easily transform 



the equation (24i to the frequency domain: 



+iV 



Xifufi) = n(/,,/ 2 ) S{fi-hlT,,f 2 -hlT 2 ) 
h,h=-N 

+N 

= UffufO J] H (fi-ji/Tuf2-ji/T z ) 



h,h=-N 



x X X(f l -i l /T l -j l /T u f 2 -i 2 /T 2 -j 2 /T 2 ) 

i] ,/ 2 =-oo 

(25) 

We know that X{f\,f 2 ) is band limited signal. So, because 



of I\(f\,f 2 ) in the equation (25i that is a lowpass filter, the 
combination of signal X(f\,f 2 ) and its shifted one, the signal 
itself is resulted. Through the explanation, equation ( |25| is 
simplified and can be written as follows: 



X(fi , f 2 ) = X(fi , / 2 )xn(/, , / 2 ) 2 ff Cfi - h ITx , h - hITi) 

(26) 

Just like the idea used in one-dimensional method, when one 
condition exists, the equation perfectly converges to x when 
N increases toward oo. 

That condition is that the equation 
muf2)ZZ 1 =-N H (fi-ji/ T iJ2-j2/T 2 ) becomes unity. 
By a little study, it is clear that this equation is true about 
all the interpolation functions, like the one that is proved 
in the one-dimensional method. Therefore, by increasing Af, 
1 goes nearer to x. But as mentioned in the last section, 
N is the number of implementation modules that because 
of practical limitations on that, it cannot be increased 
that much. In the cases with limited the amount of the 
distortion from the main signal is the difference between 1 
and n{hJ 2 )y-f h= _ N H{f,-hlT,J 2 -j 2 IT 2 ). Then, we 
can formulate the error as the equation below: 



dfidf 2 

-w, 

(27) 

In this case, the only difference between the 2-D and 1 -D 
is that £ and 23 will be matrixes and the § will be a tensor. 
With these definitions equation (21) can be written for two- 
dimensional case and it can be solved, so the compensated 
coefficients will be found. 



IV. Simulation Results and Discussion 

We utilized MATLAB® simulation environment to evaluate 
and compare the performance of methods. To have fair com- 
parison, initial band limited signals are produced randomly, 
and the performance of each method is averaged over 100 
signals. The initial signal is FFT lowpass filtered version of 
white Gaussian noise signals. To show the significance of this 
method, the sampling rate is performed at the Nyquist rate. 

The performance criterion for our simulations is the Signal 
to Noise Ratio (SNR) in dB. To avoid transient errors at the 
end points, SNR is calculated for interior points and 10% 
of the end points are ignored. As illustrated in Fig. [3] the 
SNR increases monotonically in dB for classical method as 
the number of modules increases, while the optimum method 
increases exponentially as the number of modules increases. 
This means more than 250dB for simple S&H interpolation 
and this is quite impressive in real engineering applications. 

Fig.|4]shows similar results for the Linear Interpolation (LI). 
The difference between the classical method and the optimum 
at the first few numbers of modules is not very significant. 
However, as the number of modules increases, the difference 
becomes apparent. This shows that our method can find the 
optimum coefficients for any interpolation independently. 

Classical vs Optimum (S&H) 



250 



200 



100 



50 



-A- Opt. Modular 
* - Classical Modular 



Number of Modules 

Fig. 3. SNR vs. the number of modules for Classical Modular method and 
our proposed Optimized method for S&H interpolation. 

To study the effect of noise, we added a white Gaussian 
noise to the band limited signal. This is the model of the 
electronic devices that generate thermal noise. Fig. [5] shows 
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Classical vs Optimum (S&H) for 2-D 
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Fig. 4. SNR vs. the number of modules for Classical Modular method and Fig. 6. SNR vs. the number of modules for Classical Modular method and 
our proposed Optimized Modular method for Linear interpolation our proposed Optimized method for S&H interpolation for 2-D signals. 



that for low SNR, the performance of our method is not so 
significant over classical method, however as much as the 
power of noise decrease our method give much greater SNRs. 

These simulations show that using the modular method 
by means of optimized coefficients enhance its performance 
dramatically both in noisy and noiseless environments inde- 
pendent of the type of interpolation. 

Because of the importance of the two-dimensional signals, 
similar simulations are performed on them and their results are 
shown in Fig. [6] As is clear in figure, addition of coefficients 
has boosted the speed and with fewer modules we can reach 
better results than previous one. This shows inclusiveness of 
our method in other dimensions. 

As in [ 14] and [ 15), by combination of two Modular and It- 
erative method we can reach significant results in interpolation. 
Also, here we tried to show the modular method's inclusive- 
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Fig. 5. SNR after applying methods vs. SNR of the signals before applyibg 
the Classical Modular method and our proposed Optimized method for 5 
method and the S&H interpolation. 



ness by means of compensated coefficients that are derived 
from this Hybrid Method. For this reason, we compared 
three methods, Classical Iterative (without modular), Hybrid 
Method (Iterative with Modular) and Hybrid Compensated 
Method (Iterative with Compensated Modular) with each other 
in both one-dimensional and two-dimensional domain versus 
number of iteration. As shown in figure |8(a)| and |8(b)| there 
is a significant progress in Iterative Method, so with just 
10 iteration in 1-D and 13 iteration in 2-D we can reach 
its maximum possible value and this means saving time and 
hardware in exchange for calculation of coefficients that needs 
to be calculated just once. 

Now to have a more tangible example for our method, we 
apply it to image enlargements. For this purpose we use Lena's 
known image. First after a continuous sampling the 512x512 
image is converted to 256 x 256 image. Then, this small image 
is going to be upsampled by means of various and up to 
date methods. Results are compared to the main image by 
performance criterion PSNR in dB. The results are written in 
the table [I] According to these results our method in iterative 
case has better performance than all the existing methods, 
while it has less complexity in comparison with many of these 
methods. 

For a visual comparison, the upsampled images by various 
methods are shown in figure 8. With a little precision, these 
differences are apparent in the image, specially, in the edge of 
the hat and eyes the difference is obvious. All the simulations 
show the capability of the recommended method in various 
fields. 

V. Conclusion 

A novel Optimized Modular method is proposed for com- 
pensating error of any interpolation system. We add the 
optimized coefficients calculated in a very simple manner into 
the Classical Modular method in order to maximize its perfor- 
mance as much as it could be. Also, the proposed method is 
generalized to two-dimensional domain. The simulations show 



TO BE SUBMITTED 



6 





(b) 


lEf 1 

(c) 


it \j 

(d) 







(e) (f) 
Fig. 7. Comparing various methods, (a) main image, (b) Optimized Hybrid, (c) Hybrid, (d) Iterative, (e) SAI and (f) WZP 
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Classical Iterative vs Hyprid and Optimized Hybrid (S&H) 
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TABLE I 

COMPARISON OF THE RESULTS OF LENA'S IMAGE BY VARIOUS 
METHODS VS. PSNR (DB) 



Method 


Lena 


Bilinear 


30.13 


Bicubic 


31.34 


NEDI 1 16 1 


34.10 


WZP-Haar~n6) 


31.46 


WZP-Db.9/7 |16] 


34.45 


Carey et al. 1 1 7 1 


34.48 


HMM |18| 


34.52 


HMM-SR"P9) 


34.61 


WZP-CS |16| 


34.93 


SAI f?J 


34.74 


Iterative (2 iter.) 1 1 1 1 


35.25 


Iterative (10 iter.) (tl] 


37.39 


Hybrid (2 iter, and 1 mod.) (l5l 


37.12 


Opt. Hybrid (2 iter, and 1 mod.) 


37.41 



(a) 1-D 

Classical Iterative vs Hyprid and Optimized Hybrid (S&H, 2D) 
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1 2 3 4 5 6 7 8 9 10 1 1 12 13 14 1 5 16 17 18 19 20 
Number of Iteration 



(b) 2-D 

Fig. 8. Comparison of 3 method, Classical Iterative, Hybrid Method and 
Hybrid Compensated Method vs. iteration number for S&H in (a)l-D and 
(b)2-D 

very high improvements versus classical methods. This pro- 
posed method is also more favorable in terms of computational 
complexity with respects to other error compensating methods, 
since not only the modular method has very simple algorithm 
but also we could reach better SNR values with just 2 modules 
rather than 5 modules in the classical method. In addition, 
combining this method with Iterative Method can result in a 
more powerful method and has better performance in both 
one-dimensional and two-dimensional domain in comparison 
with many other interpolation methods exist in world. 
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